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O . Abstract 

The properties of an anisotropic fluid outside a star or a black hole 
embedded in an expanding universe are investigated. One finds that, in 
Painleve-Gullstrand coordinates, the heat flux of the cosmological fluid 
vanishes, in spite of the nonzero value of a non-diagonal component of the 
stress tensor. The pressures and energy density of the fluid are regular 
at r = 2m, divergent at r = and change sign at certain values of their 
arguments. 
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1 Introduction 

Isolated black holes have been investigated in great details in the last decades. 
However, black holes (BHs) embedded in an expanding universe are also im- 
portant and even more realistic. The apparently simple question whether the 
cosmological expansion takes place even locally (at the level of an atom or of the 
Solar System) is a very complicate one and is still an unsolved problem [TJ [2] . 
Different authors gave different answers using different arguments f3J. The so- 
lution to the problem depends on the model of the universe (say, the commonly 
used Friedmann-Lemaitre- Robertson- Walker (FLRW)). 

The purpose is to combine both classes of solutions (local and cosmological) 
and to find exact solutions for the gravitational field of a compact object embed- 
ded in a cosmological background, namely a dynamical BH [H O IHl [3 [HI IH] • The 
1st author who has taken into consideration the effect of the cosmic expansion 
on local gravitational systems was McVittie 10,. He imposed a constraint on 
the (time dependent) mass of the central object for to avoid the accretion of 
the cosmic fluid into the central mass. The physical meaning of the McVittie 
solution of Einstein's equations is still under debate nowadays. A more recent 
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metric describing a spherical mass in a cosmological background is the Sultana 
- Dyer space-time . Their line element is conformal to the Schwarzschild one 
and the stress tensor corresponds to two noninteracting perfect fluids. 

Throughout the paper we use geometrical units G = c = 1, where G is 
Newton's constant and c is the velocity of light. 

2 Painleve-Gullstrand coordinates 

The Schwarzschild exact solution for the geometry outside a star or a BH 

ds^ = _(1 - —)dtl + (1 - —)-^dr^ + r^dn^, (2.1) 
r r 

revealed a series of surprises: the incompleteness of the metric covered by the 
original Schwarzschild coordinates; a highly nontrivial global structure; the dy- 
namic nature of the geometry, despite its static form revealed, for example, by 
Hawking's radiation [1^. In (2.1) ts is the Schwarzschild time and dfl^ is the 
metric on the unit 2-sphere. To get rid of the apparent singularity of the met- 
ric at the horizon r = 2m, Painleve and Gullstarnd (P-G) used the following 
temporal transformation [12l [131 HI] 

t = ts + 2V2mr + 2mln- (2.2) 

^/r + v2m 

Therefore, the line element appears as 

ds^ - -(1 - —)de + dr^ + 2^[^dtdr + r^dn\ (2.3) 
r \ r 

where t is the free-fall time, that is the proper time experienced by an observer 
who free-falls from rest at infinity. We chose the " -|-" sign in front of the square 
root in order to deal only with the inward moving free particles (along a geodesic 
curve with dr + yj2m/rdt — 0, the velocity dr/dt = — ■y/2m/r is negative). 

The geometry (2.3) is stationary, namely invariant under time translations 
(however, it is not invariant under time reversal because of the nondiagonal 
term). In addition, a constant time slice is simply flat space. We also emphasize 
that (2.3) represents physical space free falling radially into the BH at the 
Newtonian escape velocity y^2m/r. The proper time of one observer at rest 
{dr = d9 = d(t) = 0) is dr = -y/l — (2m/r)dt. Moreover, the P-G metric paints a 
picture of space falling like a river (waterfall) into a spherical BH; its 4 - velocity 
evolves by a series of infinitesimal Lorentz boosts induced by the tidal changes 
in the river speed from place to place. 



2 



3 Conformal Painleve-Gullstrand coordinates 



To look for the influence of the expanding universe on the local evolution of a 
BH, let us consider a time dependent P-G geometry given by 



(1 )dt^ +dr^ + 2 J dtdr + r^dVl? 

r \ r 



(3.1) 



where the conformal factor a{t) will be related to the scale factor in the FLRW 
universe. When r » 2m the metric (3.1) is conformally flat and therefore it 
corresponds to a spatially flat FLRW spacetime. 

Our next task will be to impose the geometry (3.1) to be a solution of 
Einstein's equations 

Gab = Rab — -^gabR'^c — ^''^Tab (3.2) 

where the Latin indices run from to 3. To reach that goal, the stress tensor 
Tab must have the nonzero components 

. 2d 4d 1 2m Sd^ 3d 1 2m 
87rT\ = — + — W , 8^r"o = T + ^\ , 

r, d^ 2d d 1 2m „ 2(2d^ — ad) 1 2m 2md 



= 8ttT% = — - — + — ^/ — , SttT^o = ^ ^ 

r \ r a-^r^ 

(3.3) 

where d{t) = da/dt > and the indices (0, 1, 2, 3) corresponds to {t, r, 9, </>). The 
expressions (3.3) have been obtained thanks to the GrTensor II package and the 
software package Maple. The scalar curvature is given by 

R\^'4-^J^ (3.4) 

It is worth noting that R\ diverges when r — >■ 0, irrespective of the form of the 
function a{t). However, the scalar curvature is regular at r = 2m due to the 
regularity of the P-G metric there. That is in contrast with the results from 
[U 121 ini m] where R\ is divergent not only at the origin r = but also at the 
event horizon. We also observe that, far from the central mass (at r >> 2m) 
all components of T1 depend only on a(t) and its time-derivative, as it should 
be. In addition, T\ = T\ = T% and T\j = 0, i.e. is diagonal. 

Let us take now a congruence of observers with the velocity vector field 

u'^=(-,--\[^,0,o], < = (3.5) 
\a a \ r I 

For an isolated BH we may put formally a{t) = 1 in (3.5) and becomes 
that is the Newtonian escape velocity. It could be checked that with 
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the ansatz (3.5) the energy-momentum tensor corresponds to an anisotropic 
fluid given by 

Tab = {P± + p)UaUb + p±gab + {Pr " P±)SaSb + Ua.qb + UbQa, (3.6) 

with p{r,t) = TabU°-u^ = —T^o the energy density of the fluid, Pr{r,t) = T\ - 
the radial pressure, p± = T\ = T\ - the pressures on the transversal directions 
9 and 0, s" is a spacelike vector orthogonal to u"-, with Sau'^ = 0, SqS° = 1, is 
the heat flux with gaw" = and it is given by the expression = —T\u^ — pu°- ^ 
obtained from (3.6). Using now (3.3), (3.5) and (3.6) we find that 

s'^ - (^0,i 0,0^ , q'^^O; -p + 2pi (3.7) 

In spite of the fact that T\ ^ 0, we have obtained a vanishing heat flux q° = 
0. That is perhaps related to the geodesic character of the congruence (3.5). 
Indeed, it is easy to calculate that = it^Vati'' = (the observers of the 
congruence are in free fall). 



4 Misner - Sharp energy 

The Misner-Sharp (MS) quasilocal energy E{t,r) [S], with its Weyl (Ew) and 
Ricci (Efi) parts, is useful to detect localized sources of gravity. In the case of 
spherical symmetry, E is obtained from [16l [171 HI] 



1 - ^ = 9^'RaRb, (4.1) 

where R = a{t)r is the areal radius and R a — dR/dxa- One finds, in the 
spacetime (3.1), that 

E(t,r) ^ malt) ^ dr^W (4.2) 

2a y r 

While the 1st term on the r.h.s. represents the Weyl energy Ew = ma{t), the 
other two terms give the Ricci energy Ef{. Keeping in mind that the cosmic 
fluid is anisotropic, Er may be written as [HI [20] 

47r 

Er = —R^p-pr+pi_), (4.3) 
With p, Pr and p± from (3.3) we find that we indeed have 

E{t,r)^^-ar'J^. (4.4) 

If the fluid were isotropic, we should have Pr = p± and, from (4.3), Eji = 
{^)R'^p, as expected. 
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Let us further calculate the total radial energy crossing a r — const, hyper- 
surface E [21] for some time interval 

W = J T^u^UaV^ dt dO d<j) (4.5) 

where n° = (0, 1, 0, 0) is the unit spacelike normal vector on E, hab = gab — naUb 
is the induced metric on E and h — det{hab)- A direct evaluation of the integrand 
gives us T'lu^Ua — 0, where the equations (3.3) and (3.5) have been used. The 
fact that = is not surprising if we remember the P-G observers are in free 
fall (the acceleration vector of' = 0) and the energy flux is vanishing. 



5 Kinematical quantities and the Raychaudhuri 
equation 

We have seen that the congruence (3.5) is geodesic so that the acceleration is 
zero. As far as the scalar expansion of the congruence is concerned, one obtains 

^ ^ n 3a 3 1 2m ,^ 

Keeping in mind that a{t) and the radial coordinate are independent variables, 
G can be made positive or negative. In addition, it is nonzero when a{t) = 
const., thanks to the fact that radial velocity is nonzero. The energy density 
and pressures from (3.3) show similar properties. This is an effect of the two 
competing factors: the stationary gravitational field of the central mass and the 
expansion of the universe. That is clear from (5.1) where the 2nd term on the 
r.h.s. is always negative and diminishes the 1st term. It is a consequence of the 
state of free fall of the ingoing fluid particles. 

For the time evolution of the expansion, (5.1) yields 

^ _ „^ Q _ — 6d^ ^ 3a ^ 3a 1 2m 9m 

2a^r V r 2a^r^ 

For an isolated BH (a = 1), we note that & — ~9m/2r^, even though the P- 
G metric is not time dependent (G < is related to our choice for the radial 
velocity, that is < 0). The shear tensor of the congruence is given by 

2ma / 2m a 1 2m 2ma 
o-QO ^ — , 0-11 = -\/ , CTio = — 5- 

(5.3) 

CT33 ar 2m 3m 

Contrary to the previous results, the components (5.3) are of the form /(a) g{r) 
and, therefore, they vanish only when r — 00 (for 0-22 and crsa, this is valid for 
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the mixed components). They become, of course, null when to = 0, because of 
the homogeneity of the space. The Raychaudhuri equation 

e + — - Vbfl'' + CJaW'"'' - UJabUj'^'' = -Rabu'^v' (5.4) 

gives the time evolution of the scalar expansion in terms of the kinematical 
parameters and the Ricci tensor of the spacetime. The vorticity tensor ujab is 
zero for our congruence and also . With 0, Q and aab from, respectively, 
(5.1), (5.2) and (5.3), we see that (5.4) is satisfied, with RabU°'u^ obtained from 

Rabu'^u'' ^8n{p+^Tl). (5.5) 

We observe also that all the kinematical quantities are divergent at the singu- 
larity r = 0, irrespective of the value of a{t). 



6 de Sitter embedding of the P-G metric 

Once we have analyzed the time dependent P-G space, we intend now to choose 
a particular expression of the scale factor a{t) and study that case in detail. 
We shall follow the same procedure as in [15]: take firstly a{t) = e^* and then 
change to = —{1/H)e^^^ (i > 0, — < < 0). We therefore embed the 
P-G black hole in a time dependent de Sitter spacetime and the metric is 



- 1 )dT]^ + dr^ + 2\ dr]dr + r^dn^ 

r \ r 



(6.1) 



with a{'i]) = 1/H\ri\. For r >> 2m, (6.1) acquires the conformally fiat de Sitter 
form. In contrast, with 0(77) = 1 one obtains the stationary P-G line element. 
With the above special value for the scale factor 0(77), the energy density and 
pressures of the anisotropic fiuid obtained from (3.3) are given by 




Pr = 



PA- 



3i/2 / 




/¥ 
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Stt \ 


^ 3r \ 




3ff2 
Svr 




/¥ 





(6.2) 



with p+Pr —p± — 3H'^ /Stt, that may work as an equation of state. We note that 
asymptotically (r >> 2to) the fiuid becomes isotropic {p^ = p±) and p = —pr, 
as expected for a de Sitter spacetime. The same effect is obtained if we put 
77 = in (6.2) (namely, t — >■ 00). Moreover, p, pr and p± remain finite at 
r — 2m. We also observe they change sign at certain value of r and constant 77 
or at some value of 77 and constant r. In other words, the weak energy condition 
is not satisfied for all (77, r) pairs. 
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Similar properties we find for the expansion Q and the scalar curvature 

While Q and i?"^ become constant at r >> 2m or 77 0, SttTq = 2mH^\ri\/r^ 
is negligible. One should be so for to recover the pure de Sitter values where & 
and R\ are indeed constant {3H, respectively 12H^) and the stress tensor is 
diagonal. 

7 Conclusions 

We investigated in this paper a compact object (a star or a black hole) in 
Painleve-GuUstrand coordinates and immersed in an expanding universe, pay- 
ing special attention to the influence of cosmology on the local physics. The 
anisotropic fluid is freely falling at Newtonian escape velocity and the chosen 
congruence of observers is geodesic. Therefore, there is no heat flux despite the 
nonzero value of the T q component of the energy momentum tensor. Excepting 
the shear tensor, the other parameters of the fluid (energy density, pressures, 
scalar expansion) change their sign at some value of r at constant rj, or of 77 at 
constant r, a property that is difficult to have a reasonable physical interpreta- 
tion (see, for example where the authors reached similar conclusions). 

References 

[1] C.-Y. Sun, Commun. theor. Phys. 55, 597 (2011) (ArXiv: 1004.1760 [gr-qc]. 
[2] S. A. Klioner and M. H. Soffel, ArXiv: astro-ph/04l"T363 
[3] W. B. Bonnor, Gen. Relat. Grav. 32, 1005 (2000). 

[4] V. Faraoni and A. Jacques, Phys. Rev. D76, 063510 (2007) (ArXiv: 
0707.1350 [gr-qc]). 

[5] M. Carrera and D. Giulini, Phys. Rev. D81: 043521 (2010) (ArXiv: 
0908.3101 [gr-qc]). 

[6] H. Saida, T. Harada and H. Maeda, Class. Quant. Grav. 24: 4711 (2007) 
(ArXiv: 0705.4012 [gr-qc]). 

[7] M. Anderson, J. Phys., Conf. Ser. 283, 012001 (2011). 

[8] V. Faraoni, A. Z. Moreno and R. Nandra, ArXiv: 1202.0719. 

[9] V. Faraoni et al, Phys. Lett. B671, 7 (2009). 

[10] G. C. McVittie, Mon. Not. R. Astr. Soc. 93, 325 (1933). 



7 



[11] J. Sultana and C. C. Dyer, Gen. Relat. Grav. 37, 1349 (2005). 

[12] P. Kraus and F. Wilczek, ArXiv: |gr-qc/9406042 

[13] G. Tian, S. Wang and Z. Zheng, ArXiv: | gr-qc/0509030| 

[14] E. Poisson, A relativist's toolkit "The mathematics of black hole mechan- 
ics", Cambridge University Press (2004). 

[15] H. Culetu, ArXiv: 1201.3769 [gr-qc]. 

[16] A. Nielsen and D.-h. Yeom, Int. J. Mod. Phys. A24: 5261 (2009) (ArXiv: 
0804.4435 [gr-qc]). 

[17] M. Carrera and D. Giulini, ArXiv: 0810.2712 [gr-qc]. 

[18] L. Herrera et al.. Int. J. Mod. Phys. D18: 129 (2009) (ArXiv: 0804.3584 
[gr-qc]). 

[19] B. V. Ivanov, ArXiv: 1103.4225 [gr-qc]. 
[20] H. Culetu, ArXiv: 1202.4296 [gr-qc]. 

[21] D. C. Guariento and J. E. Horvath, Gen. Relat. Grav. 44, 985 (2012). 
[22] H. Culetu, ArXiv: 0711.0062 [hcp-th]. 



8 



